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ABSTRACT 

Galaxy distances and derived radial peculiar velocity catalogs constitute valuable datasets to 
study the dynamics of the Local Universe. However, such catalogs suffer from biases whose effects 
increase with the distance. Malmquist biases and lognormal error distribution affect the catalogs. 
Velocity fields of the Local Universe reconstructed with these catalogs present a spurious overall 
infall onto the Local Volume if they are not corrected for biases. Such an infall is observed in the 
reconstructed velocity field obtained when applying the Bayesian Wiener-Filter technique to the raw 
second radial peculiar velocity catalog of the Cosmicflows project. In this paper, an iterative method 
to reduce spurious non-Gaussianities in the radial peculiar velocity distribution, to retroactively de¬ 
rive overall better distance estimates resulting in a minimization of the effects of biases, is presented. 
This method is tested with mock catalogs. To control the cosmic variance, mocks are built out of 
different cosmological constrained simulations which resemble the Local Universe. To realistically 
reproduce the effects of biases, the mocks are constructed to be look-alikes of the second data release 
of the Cosmicflows project, with respect to the size, distribution of data and distribution of errors. 
Using a suite of mock catalogs, the outcome of the correction is verified to be affected neither by the 
added error realization, nor by the datapoint selection, nor by the constrained simulation. Results are 
similar for the different tested mocks. After correction, the general infall is satisfactorily suppressed. 
The method allows us to obtained catalogs which together with the Wiener-Filter technique give 
reconstructions approximating non biased velocity fields at 100-150 km s ' (2-3 h *Mpc in terms of 
linear displacement), the linear theory threshold. 

Key words: methods; numerical, techniques: radial velocities - large-scale structure of universe 


INTRODUCTION 


lodes, sheets, filaments and voids compose the Large Scale Structure 
if the Local Universe. Understanding the apparent complexity of grav- 
tational motions that results in the formation of such an environment 
onstitutes a real challenge. For instance, progress has been made over 
he past few years to better understand the mo tion of our Galaxy in this 
ontext ( ITu11vI|200^ I 2 OO 8 I : |Tu11v et al.ll2008l) . Yet, in an era of ”preci- 
ion cosmology”, only 30% of the motion of the Milky Way relative 
3 the Cosmic Microwave Background is understood exactly. Catalogs 
of accurate distances and radial peculiar velocities of galaxies consti¬ 
tute a valuable source of information to study our neighborhood. Pe¬ 
culiar velocities are direct tracers of the underlying gravitational fields 
as they account for the total amount of matter. This is an advantage 
over redshift surveys which, although they are more easily obtained, 
acco unt solely for the luminous matter and thus are subject to biases 
(e.g. lBau^l2013h . For that reason, Cosmicflowfl whose goal is to 
map the Local Universe, develops such radial peculiar velocity cat¬ 
alogs. The final aim is to reconstruct overdensity and velocity fields 
of our neighborhood using exclusively peculiar velocities which are 
highly linear, correlated on large scales and not affected by the lumi- 
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nous bias. A large amount of observational data is necessary to ob¬ 
tain accurat e galaxy dis t ances from which the radial peculiar velocities 
are derived. ITuIIv et aP ( 1201 3h released the second generation catalog, 
cosmicflows-2, of the project which reaches 150 h^'Mpc and contained 
more than 8,000 galaxy distance and peculiar velocity measurements. 
Although superior in size and extent to the first catalog of the project, 
this catalog is not, however, immune to biases when considered without 
special treatments as stipulated by Tully et al. 2O1!;0. At such a depth, 
errors can be larger than peculiar velocities, overwhelming their signal. 
When reconstructing the full 3D velocity field, applying the Bayesian 
Wiener-Filter technique (WF, IZaroubi et al.ll 19951 1 19991) assuming by 
default a symmetric error distribution, the ne gative tail (asym metric 
with respect to the positive tail, highlighted by ItuIIv et ^ .l2013h in the 
distribution of radial peculiar velocities and Malmquist Biases give rise 
to a major spurious infall onto the Local Volume. This paper presents 
a detailed analysis of these biases and, in addition, proposes a method 
to retroactively derive overall better distance estimates by limiting the 
asymmetrical bias in radial peculiar velocity catalogs. This method re- 


^ ”[...] although individual distances may be unbiased, the homogeneous and 
non-homogeneous Malmquist effects and error bias can generate spurious arti¬ 
facts in vetocity fietds. These issues must be addressed if velocities are used to 
infer the distribution of matter.” 
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suits in the suppression of the large infall and thus contributes to mini¬ 
mize the biases responsible for this spurious general motion. 

In Section 2, the biases are explained in detail. In Section 3, the 
cosmicflows-2 features - in terms of size, data distribution and error 
distribution - are analyzed in order to mimic them with mock catalogs 
that reproduce as realistically as possible the biases. In Section 4, mock 
catalogs similar to the observational dataset are built out of reference 
cosmological simulations which resemble our Local Universe to limit 
the cosmic variance. A correction method is derived and tested on these 
mock catalogs. Three types of mocks are built to check the method 
outcome against the added error realization, the data selection and the 
constrained simulation. Namely, first type mocks differ only by the ran¬ 
domly added error, while mocks of the second type are constituted of a 
different datapoint selection from the same reference simulation. As for 
third type mocks, they are built out of different constrained reference 
simulations. After correction these mocks are input into the Wiener- 
Filter technique to reconstruct overdensity and velocity fields of the 
reference constrained simulations. Resulting reconstructions are com¬ 
pared with reconstructions obtained with their biased (mocks with er¬ 
rors and without correction) and original (mocks without errors) coun¬ 
terparts, and the effectiveness of the method is evaluated. Last but not 
least, the method is applied to the raw observational, cosmicflows-2, 
catalog. 

Throughout this paper, distances are in h^'Mpc. All figures 
are presented after a Gaussian smoothing of (over)density fields at 
2 h“‘Mpc, the linear theory intrinsic floor value, as the WF - linear 
minimal variance estimator - reconstruction technique used in this pa¬ 
per is linear by definition. Conversions between velocity fields and 
displacement fields ( km s“' t o h~‘Mpc) are ob tained using the lin¬ 
ear Zel’dovich approximation llZel’dovichlll970h and the definition of 
the peculiar velocity (see the Appendix for a detailed equation). As 
constrained simulations were obtai ned within the seven- year Wilkin¬ 
son Microwave Anisotropy Probe llKomatsu et alJ I 20 I 1 I 1 framework, 
reconstruction tests are conducted within the same framework. 


2 DEEPER CATALOGS: STRONGER BIASES ’ EEEECTS 


Distance estimates are subject to severe systematic biases which af¬ 
fect the derived radial peculiar velocities. Whilst generally all gathered 
under the term "Malmquist Bias”, three types of Malmquist Bias can 
in fact be distinguished. In addition to these biases, there is a lognor¬ 
mal error distribution which requires some attention. It is necessary to 
add that the biases are mentioned and explained considering the obser- 
va tional radia l pecul iar velocities catalog, cosmicflows-2, as presented 
in iTullv et alJ lIZOLlh . That is to say that different techn iques over the 
years, such as those detailed in IStrauss & Willi^(ll99.5h relying on un¬ 
derlying models to correct distances as a whole before deriving radial 
peculiar velocities, constitute other options. Still radial peculiar veloc¬ 
ities in the cosmicflows-2 catalogs were derived without considering 
all the biases as a choice to privilege individu al distance measur ements 
over an ensemble of distance measurements dTullv et al.ll20L3h . In the 
four following points, the different biases are disentangled to identify 
those left in the cosmicflows-2 observational catalog which lead to a 
spurious infall onto the reconstructed Local Volume. 


• The most often mentioned Malmquist Bias, referred to as Problem 
I, Selection Effect/Bias, ”r against V”, Distance-dependent, Frequen- 
tist. Calibration pr oblem, M-bias of the second kind (Kaptney , 1914; 
Malmquist. 1922: lHanlll9p: ISandagelll994 lTeerikorDilll99l [19^ 


Il99d: iHendrv & Simmonslll994 IWillicklll994 lTeerikorpilll995ir 7s 
analogous to a selection effect in magnitude resulting in underesti¬ 
mated distances. A magnitude limit in a selected sample results in a 
mean apparent magnitude smaller than it should be. As the limit de¬ 
creases, in unit of magnitudes with the distance, the bias increases. 


Namely, dwarfs and other dimmer galaxies are more and more under¬ 
represented in the sample with respect to brighter galaxies as the ob¬ 
server looks farther and farther. Consequently, their contribution to the 
mean apparent magnitude is increasingly reduced. Thus, the mean in¬ 
creases slower than it would have, with the distance or redshift, would 
all the galaxies have really been included. Such a bias is usually already 
taken care of when calibrating distance indicat ors. This is the case with 
the observational cosmic flows-2 catalog (e.e. ITuIIv & Courtoisll2012l ; 
ISorce et"^ .l20Lll2014bh . 

• The second bias is called Homogeneous Malmquist Bias and 
gathers terminology such as Problem II, General Malmquist Bias, 
Geometry Bias, ”V against r”. Classical, Bayesian, Inferred-distance 


Lvnden-Bell et al 

19881; Haul 1992; Teerikoroil 1994 Sandaed 

1994 

Teerikomi 19931 

I 990 I: iHendrv & SimmonsI 

1994 TeerikoiDil 

1995 

Strauss & Willick 

|l99.5h. Because our observations are restricted to a 


sphere centered on us, from the homogeneity of the Universe, the num¬ 
ber of observable galaxies increases with the distance. Then consider¬ 
ing a true distance (r), there is a higher probability to observe a galaxy 
at (r) -r dr than at (r) - dr. Consequently, the probability for a galaxy 
which distance estimate is r to have been located closer is higher than 
the opposite. Namely, distances are more likely to be underestimated. 
In practice, it is even more complicated than that because of the pres¬ 
ence of small scale structures. The bias is additionally a function of the 
direction of observations leading to the third Malmquist Bias, 

• The last of the Malmquist Bi as is referred to as the In¬ 
homogeneous Malmq uist Bias (e.g. iDekell 1 1994 iHudsonI 1 1994 
iLandv & Szalaviri992h . It is the result of the small scale structures, 
namely of the fluctuations of galaxy numbers. To illustrate this bias, let 
us consider a group of galaxies at (r) with null radial peculiar velocities 
for simplicity. Because of observational uncertainties, these galaxies 
are randomly scattered to the foreground and background of (r). For 
all galaxies at the same redshift, their estimated radial peculiar veloci¬ 
ties obtained on either side of (r) result in an inaccurate infall towards 
(r) giving birth to spurious structures and flows. In other words, galax¬ 
ies are more likely to be scattered from high density regions towards 
low density regions than the opposite. 

• In addition to the Inhomogeneous Malmquist Bias, the fractional 
errors on distances have a lognormal distribution. The distribution of 
fractional errors is not symmetric when considering two galaxies one 
located closer to the observer than it should be and the other one po¬ 
sitioned farther away. This asymmetry introduces a strong bias in the 
final recovered velocity field. In other words, distances, d, of galax¬ 
ies are derived from moduli (ju, apparent magnitude m minus absolute 
magnitude M) themselves obtained with observational data. Thus, as¬ 
suming that errors in magnitudes are distributed evenly for moduli, 
since a logarithmic function is used to derive distances in megaparsecs 
from moduli yu in magnitudes, errors in megaparsecs are not distributed 
symmetrically around distance values, d as shown by equations [T] 


( 1 ) 


ji = m — M 

Afi = VA)7j- -I- AM^ sj AM symmetric around /j 

d = 10<^-25)/5 

Ad = (Ayu X In(lO) X d)/5 asymmetric aroundd 

Since radial peculiar velocities, Vpec and their uncertainties are de¬ 
rived from distances and their errors by the classical equations |2l the 
asymmetric distribution propagates onto peculiar velocities. A posteri¬ 
ori, looking for a cause of the observed negative tail in the distribution 
of radial peculiar velocities, the assumption of a Gaussian distribution 
of errors on distance moduli is validated. 


^pec ^mod Hq X d 

Avpec = HoX Ad 


( 2 ) 
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where is the velocity with respect to t he Cosmic Microw ave Back¬ 
ground corrected for cosmological effects jTullv et al.ll2013h . 

This asymmetrical error distrihution results in a bias in peculiar ve¬ 
locities. Large negative peculiar velocities have an abnormal large error 
which cannot be properly propagated in the WF because the sign of the 
error on a particular distance modulus is unknown. 

As an example, let us consider a galaxy at 100 Mpc, with 
Vpec=0 km s“*, Ho=75 km Mpc“* s“'. Assuming that the distance 
modulus of this galaxy is obtained with the Tully-Fisher relation (the 
main provider of distances in cos micflows-2), the 1-cr uncertainty on 
the measurem ent is about 0.4 mag jTullv & Courtoisl2012l : [^rce et alj 
l20lll2014bh . A 1-(T error on the distance modulus gives either 34.2 or 
35.8 mag instead of yu = 35 mag. The derived distance is either 69 or 
145 Mpc, giving a radial peculiar velocity of 2311 or -3338 km s“*. As 
a result, a larger error is made by allocating a radial peculiar velocity 
of -3338 km s^' to this galaxy than when assuming 2311 km s“'. Flow- 
ever, the sign of the error on the distance modulus is unknown, hence, 
although -3338 km s“* is assigned to this galaxy, proportionally the 
same uncertainty is attributed to that peculiar velocity as it would have 
been if the value 2311 km s“‘ was given to that same galaxy. Since, the 
WF technique uses errors on peculiar velocities as an indication of the 
signal strength in the auto-correlation matrix, these velocities would 
be proportionally identically weighted despite the absolute error of the 
negative value being greater than the positive one. This asymmetrical 
bias contributes, together with the uncorrected for Malmquist biases, 
to the spurious overall large infall seen in the reconstructed velocity 
field. 


Reducing such biases is a complicated task. Grouping galaxies in 
groups and clusters m ight reduce these bia s es as recommended by e.g. 
iBurstein et alJ jlQQOh and iNusser & Davlsl j201 ill : because galaxies in 
the same cluster or group approximately share the same distance, aver¬ 
aging over galaxy distance estimates (to keep only distance estimates 
of groups and clusters) reduces errors on distances, thus on radial pecu¬ 
liar velocities, of clusters and groups by the square root of the number 
of estimates. Retracin g the footsteps of methods prop osed for instanc e 
bv IPekel et alJ jl999h for the POTENT technique or bv IlTudsor] jl994ti . 
in addition to grouping, distances can be corrected or more appropri¬ 
ately, their uncertainties can be minimized. A method to reduce un¬ 
certainties on distances is presented in the following sections. Unlike 
the above mentioned methods, which assumed an u nderlying distance 
distri bution as a prior to correct distances (see also lLvnden-Bell et alJ 
Il988h . this paper proposes that a velocity distribution serves as a basis 
for the correction. Both the lack of knowledge regarding the true under¬ 
lying distance distribution and the multiple distance-indicator nature of 
cosmicflows-2 lead us to this choice. The method is thus based on cor¬ 
recting peculiar velocities first to decrease the flatness and skewness in 
the radial peculiar velocity distribution. The resulting distances derived 
retroactively from these corrected radial peculiar velocities are by ex¬ 
tension better distance estimates. The set of corrected radial peculiar 
velocities and corrected distances results in reconstructions where the 
spurious infall due to biases is drastically decreased. 


3 ANALYSIS OF COSMICFLOWS-2 CATALOG 


Cosmicflows-2 is the second generation catalog of galax y distances 
built by the Cosmicflows collaboration. Published in ITuIIv et all 
j2013h . it contains more than 8,000 accurate galaxy peculiar ve¬ 
locities. Distance measuremen ts come mostly from the Tully-Fisher 
relation j Tullv & Fisheil fl977li a nd the Fundamental P lane methods 
jColless et al.l 200 ih. Clepheid s ( Freedman et al.l 1200 ih . Tip of the 
Red Giant Branch jLee et alJll993h , Surfa ce Brightness F luctuation 
jTonrv et alJl200lh . supemovae of type la j.Tha et alJl2007h and other 
miscellaneous methods also contribute to this large dataset though to 
a minor extent (~ 12%). The Cosmicflows team put lots of effort into 


gathering these data and ensuring that they are on a common scale, 
proceeding from closer to farther scales to check the consistency of the 
data and adjusting zero points when necessary. They are able to do so 
thanks to the numerous measurement overlaps between the different 
techniques. Finall y, they polish their scale using that obtained with the 
Spitzer Telescope jSorce et ^l2012h which gives the great advantage, 
among others, of a 1% all sky consistency erasing any concern of dis¬ 
crepancies between the northern and the southern observations. From 
these considerations, cosmicflows-2 is assumed hereafter to consist of 
distance measurements on a common scale. 

As noted in the previous section, grouping may reduced biases. 
Therefore, in this paper, the last grouped version of cosmicflows-2 
(Tully, private communication) is used as a reference to reproduce as 
realistic as possible biases effects in a radial peculiar velocity catalog 
of our neighborhood. 552 groups and 4303 single galaxies can be iden¬ 
tified in the dataset reducing the number of radial peculiar velocities to 
4855. This number is stil l large enough f or the c urrent scientific pur¬ 
pose as demonstrated by IPoumler et alJ j2013bll3 9i. Figure [T] shows 
the spatial distribution of these groups and single galaxies in a ± 5 
h^'Mpc thick slice of the supergalactic XY plane on to p of galaxies 
from the 2MASS redshift catalog within the same slice jHuchra et alJ 
l2012h . This distribution is inhomogeneous at least on two important 
respects: 1) the density of measurements decreases with the distance to 
the observer (for instance, 98 % of the data are within a 160 h^'Mpc ra¬ 
dius sphere centered on us but half of them are within a much smaller 
volume, a 60 h^'Mpc radius sphere centered on us), 2) a large zone 
devoid of data due to galactic extinction is visible. In order to mimic 
as much as possible cosmicflows-2, and as a result biases, in the re¬ 
lated validation tests, these two features must be well reproduced by 
the mock catalogs. A feature which should also be mimicked comes 
from the error distribution. This distribution is strongly bimodal with 
a peak around 8-10% fracti onal uncertainty o n distances and another 
one at approximately 18 % jTullv et al.ll20L3h . This bimodality is due 
to the different distance indicators and to the grouping: measurements 
obtained with supemovae, tip of the red giant branch and cepheids are 
given smaller uncertainties than measures obtained with Tully-Fisher 
or Fundamental Plane relations. As for grouping, averaging over a high 
number of measurements reduced the uncertainty on the measure of the 
group with respect to the estimate for a single galaxy. 


4 DERIVING A CORRECTION METHOD, TESTS ON 
MOCK CATALOGS 

Instead of correcting distances, as widely proposed in previous meth¬ 
ods (e.g. lPekel et alJ[l999l : lHudsonlll994h . this paper proposes first to 
adjust peculiar velocities and then go back to correct the distances. 
Consequently, the process is based on the distribution of radial pecu- 
liar velocities rather than on the radial distribution of galax ies (e.g. 
iLandv & Szalavll992HHudsonll994h . lSheth & Diaferi^ l 200 ih proved 
that the distribution of radial peculiar velocities considering groups and 
clusters (namely removing virial motions) should be a Gaussian. Con¬ 
sequently, unless the Milky Way is at a peculiar position in the Uni¬ 
verse, which is highly improbable, the distribution of radial peculiar 
velocities obtained from our position should be close to a Gaussian 
too. As a matter of fact, Gaussianity is found in mock peculiar ve¬ 
locity catalogs drawn from N-body simulations. Namely, dark matter 
halos (accessible with dark matter only simulations, themselves well 
described by the linear theory) are equivalent to groups, clusters or iso¬ 
lated galaxies. To derive a method to minimize biases in observational 
datase ts, this Gaussianity will be the initial assumption. iBhattacharval 
j2008h noted, when studying radial peculiar velocities in simulations, 
that although the overall distribution of radial peculiar velocities was 
confirmed as Gaussian, we still can be affected by: 
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Figure 1. Distribution of groups and single galaxies, in a ± 5 h“*Mpc thick 
slice, in cosmicflows-2, observational catalog of the Cosmicflows project (green 
dots) superimposed to the 2MASS redshift catalog (tiny black dots, same slice). 
This view (supergalactic XY plane) shows clearly the Zone Of Avoidance. 


• the cosmic variance due to the particularity of our neighborhood, 

• Poisson noise due the restricted size of the sample of peculiar 
velocities. 

Still, a major advantage of this paper study comes from the fact that: 

• constrained simulations of the Local Universe are used to produce 
a set of mocks on which the method is tested, 

• mock catalogs mimic as much as possible the characteristics of 
the observational catalog under study. 

As a result, both cosmic variance and Poisson noise are reduced. In 
any case, peculiarities of the Local Universe make purely mathematical 
statistical methods not optimal. 


4.1 Mock catalogs of peculiar velocities 


To limit the impact of the cosmic variance in the following tests, sim¬ 
ulations resembling the Local Universe down to the linear thresh¬ 
old (2 h“‘Mpc) are chosen. Our simulations have 512^ particles 
and are 320 h^'Mpc wide. They were computed within the frame¬ 
work of the CLUE5 1 proj ect (Constrained Local Unive rse Simula¬ 


tions, 


Gottlober et al 


201011 using the method described in ISorce et alj 
1 I 2 OI 4 X Unlike typical simulations, constrained simulations stem from 
a set of constraints that can be either redshift surveys or radial peculiar 
velocities. Simulations used here obey a set of radial peculiar veloc¬ 
ities. Most importantly, a look alike for all the major structures and 
voids of the Local Universe can be found in these simulations ensuring 
that the Large Scale Environment is similar to our neighborhood. Sev¬ 
eral realizations are tested to measure the robustness and the accuracy 
of the method discussed in the next subsection. Results are similar, 
if not identical, for every mock tested. Consequently, while the hgures 
are found by averaging over the entire mock set, one mock (chosen ran¬ 
domly) built out of one of the aforementioned simulations is presented 
in this paper in terms of plots. Considering the box to be equivalent 


^ http://www.clues-project.org/ 
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Figure 2. Example of a distribution of selected halos (blue dots) in the XY 
plane (± 5 h“*Mpc slice in Z) of a cosmological simulation which resembles 
the Local Universe. Halos are selected to build a mock similar to the second 
catalog of the Cosmicflows project. Black and green contours show the density 
field and the mean density of the cosmological simulation. 


to the Local Universe, we place an observer at the center of the box, 
coordinates can be defined similarly to observational supergalactic co¬ 
ordinates. The XY plane in this set of coordinates, is shown in Fig¬ 
ure |2] a Shapley supercluster candidate in the top left comer, a look- 
alike of Coma in the top middle, and a candidate for Perseus-Pisces 
around [-25,25] h^'Mpc, Virgo’s look-alike is close to the center and 
the Centaurus-Gre at Attractor region-like is on Virgo’s left side. Using 
Amiga halo finder jKnollmann & Knebell2009t) . a list of halos is drawn 
from this simulation. Halos are then selected and prepared to match 
cosmicflows-2 grouped catalog in three steps: 

• We seek to have a similar repartition of data points (number, spa¬ 
tial coverage and distribution including the Zone Of Avoidance). Sub¬ 
sequently, every halo in a zone similar to the Zone of Avoidance is 
removed from the list. This zone is defined as a cone with the apex at 
the center of the box (where the observer is assumed to lie) and assum¬ 
ing the same orientation within the XYZ volume as the observational 
one in the supergalactic XYZ volume. Then, an histogram with bin 
size of 20 h^'Mpc is derived for the observational cosmicflows-2 cata¬ 
log providing the number of measurements in each bin or 20 h^'Mpc 
slice. For each 20 h^'Mpc slice, the same number of halos as found 
in cosmicflows-2 is selected randomly in the list of halos. In case of 
multiple halos in the same region, the most massive halo is privileged 
to mimic the grouping applied to cosmicflows-2. On Figure |2l the re¬ 
sulting compiled list of halos is visible as blue dots in a 10 h~'Mp c 
thick slice in the XY plane. As predicted bv isheth & Diaferiol (1200 ih . 
before inserting errors, the halo radial peculiar velocity distribution 
(computed with respect to the box center, where the Milky-Way like 
is assumed to be) of this mock catalog can be modeled by a Gaussian 
visible in Figurel^in blue. This mock is called hereafter original in the 
sense that it has no error on galaxy distances yet. 

• Next, we add errors to radial peculiar velocity distributions pre¬ 
serving the asymmetry problem. Accordingly, a Gaussian distribution 
of errors with 0.17 magnitude scatter (as on average by assumption in 
the observational catalog) is added to distance moduli (not to distances) 
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Figure 3. Distribution of radial peculiar velocities in the original (without en'or) 
mock catalog (blue solid thick histogram) and in the biased mock (black dotted 
histogram). The original distribution of radial peculiar velocities can be mod¬ 
eled by a normal distribution (blue solid light curve). The distributions of cor¬ 
rected and 1.5 corrected radial peculiar velocities are shown by the red dashed 
thick and green dot-dashed thick histograms. Gaussians can also fit these distri¬ 
butions (red dashed light and green dot-dashed light lines). 


Figure 4. Radial peculiai' velocity distributions in CF2-Biased (black dotted 
line), CF2-Con'ected (red dashed thick line) and CF2-1.5-Con'ected (green dot- 
dashed thick line) catalogs. A larger tail is visible in the negative side of this 
diagram due to the error bias in the velocity distribution of CF2-biased. After 
con'ection, the distributions of radial peculiar velocities in cosmicflows-2 cor¬ 
rected datasets can be fitted by Gaussians (red dashed and green dot-dashed 
light lines) with parameters close to the theoretical ones (obtained with mocks). 


and distorted distances and corresponding radial peculiar velocities are 
computed. Because distances are distorted the Malmquist biases are 
also reproduced. Figure 0 displays the distribution of radial peculiar 
velocities for this mock with errors, hereafter biased mock, by a black 
dotted line. This distribution is flatter than a theoretical Gaussian with 
a slightly larger tail on the negative side confirmed by a negative skew¬ 
ness value. This distribution and that of cosmicflows-2 represented also 
by a dotted black line but on Figure[4]are similar to each other confirm¬ 
ing the 0.17 mag chosen scatter. 

• Finally, the inserted errors are the ’true errors’ but are unknown to 
observers. Observers have at their disposal only the l-cr uncertainties 
which differ on the distance indicators and the grouping. This frac¬ 
tional distributi on of l-cr uncert ainties is globally highly bimodal in 
cosmicflows-2 llTullv et alJl^l3h . accordingly distances are assigned 
in majority a 17% and in minority a 5% l-cr fractional uncertainty. 
This assignment is not done entirely randomly but in accordance with 
the real inserted errors in the previous step, analogously to how ob¬ 
servers attribute the fractional uncertainty in function of the distance 
indicator. 

We repeat this three-step process several times with this and all 
the other constrained simulations at our disposal. F-tests show that re¬ 
sulting radial peculiar velocity distributions are always similar at the 
98 ± 2 % confidence level. 


4.2 Method: minimization of the Biases using Mocks 

We intend to present a method to correct overall distance estimates by 
correcting first radial peculiar velocities. Figure |4] reveals a skewness 
towards negative peculiar velocities in cosmicflows-2 confirm ing the 
asymmetry phenomenon. Flowever, ISheth & Diaferiol bOOlh proved 
that the radial peculiar velocity distribution should be a Gaussian. Con¬ 
sequently, to derive a correction for observational datasets, this Gaus- 
sianity is the initial assumption in the method proposed in this paper. 

After application of the correction-process, catalogs should have 
a Gaussian radial peculiar velocity distribution with variance and loca¬ 
tion peak similar to that of the original ones. In mock cases, variance 


and peak are determined by mocks with original positions, hence ve¬ 
locities. Regardless, in the different original mocks, standard deviations 
of Gaussians fitted to radial peculiar velocity distributions are always 
around 300 ± 50 km s“'. Actually, F-test results, given on average 
in Table [T] column 6, reveal that a Gaussian of variance 300 km s“* 
can be accepted on average as a model for original radial peculiar ve¬ 
locity distributions at the 98 ± 2 % confidence level. As a result, a 
biased peculiar velocity is modified according to its probability of be¬ 
longing to the theoretical Gaussian (with a typical standard deviation 
of 300 km s“*) and according to its uncertainty. Two cases can be 
distinguished, either the radial peculiar velocity is positive or it is neg¬ 
ative. Then corrected radial peculiar velocities {vpec c) are derived with 
equations 0 and |4] devised in this work. 

if Vp,, > 0, Vp,,, = w [pivp,, - A) -I- (1 - p)ivp,, A)] + (1 - ui) Vp,, (3) 

if Vp^i- < 0, Vpgf-f- — w [pfHyjec (1 ~ P)iPpec ~ ^)] (1 ~ Vpgf- (4) 

where A is the radial peculiar velocity uncertainty {ISvpec)- p and w are 
two weight factors between two quantities. More precisely, p is the 
probability that a radial peculiar velocity does not belong to the theo¬ 
retical Gaussian (thus it needs to be corrected and it should either be 
reduced if it is highly positive or increased it is highly negative), p 
acts as the weight factor between Vp^^ — A and Vp^c -l- A to establish the 
quantity by which a velocity is decreased/increased. This gives an in¬ 
termediate velocity value but no information regarding the uncertainty 
on the velocity estimate has yet been considered. To use this additional 
information, w is the weighted uncertainty in opposition with A (’un¬ 
weighted’ uncertainty), w is the weight between the intermediate value 
and the measured value. The weight factors are derived with the fol¬ 
lowing equations: 

n-ll>pecl pi ^ A 

p = 1- I expi - -)dv w = —x(max( —)-t0.08) (5) 

J-oo Hod Hod 

The 0.08 term in w is required to ensure that we keep a minimum of 
trust even towards radial peculiar velocities with the maximum (max) 
fractional uncertainty. Several values between 0 and 1 were tested, 0.08 
being found to be the best parameter to retrieve a distribution close to 
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Figure 5. Histograms of fractional errors on distances (top) and absolute enors 
on radial peculiar velocities (bottom) in the biased mock (dotted black lines) and 
in the corrected mock (solid red lines). Absolute fractional errors on corrected 
distances have a median at 0.03 and a mean at 0.04. 


the theoretical Gaussian (at the 80 ± 3 % level) at the end of the cor¬ 
rection process. In other words, to correct velocities, their uncertainty 
(weighted, w or not, A) and their ‘position’ on the theoretical Gaussian 
{p) are the two parameters used. 

As A corresponds to the 1-cr uncertainty, an iterative process 
was necessary. Namely, after the first iteration, peculiar velocities 
which are still in the l-cr and beyond range of the theoretical Gaussian 
distribution are input again in the equations. After, this second 
iteration, velocities in the 2-a and beyond range are corrected again. 
Four iterations are sufficient to retrieve the Gaussian distribution 
without any long (negative and positive) tails and skewness. For 
consistency, corresponding distances need to be computed to correct 
for the Malmquist biases. They are derived using corrected radial 
peculiar velocities and the classical formula|2]in reverse order. 

After correction, the distributions of radial peculiar velocities, 
fractional errors on distances, and absolute errors on radial peculiar 
velocities are all shown on Figures[3]and|5]in red for the mock selected 
for display. Figure also presents, before correction, the histograms 
of fractional errors on distances (top) and of absolute errors on radial 
peculiar velocities (bottom) by black dotted lines. Absolute fractional 
errors on corrected distances have a median at 0.03 and a mean at 0.04. 
Fractional errors on distances are distributed on an approximate Gaus¬ 
sian and the distribution of uncertainties on radial peculiar velocities 
is less fiat and contained in a smaller interval of values. These distri¬ 
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Figure 6. Top: Residual between biased (black), corrected (red) distances and 
true distances versus true distances. The standard deviation is about 15 h“*Mpc 
before correction while it is only about 4 h“*Mpc after con'ection. Bottom: Bi¬ 
ased (black) and con'ected (red) velocities versus true velocities. The thin black 
dashed line shows the ideal correlation. The thick black dash-dotted line is the 
correlation before correction while red lines stand for correlations after correc¬ 
tion. The thinner, triple dots-dashed, red line shows the correlation without as¬ 
suming any special errors on velocities while the thicker dotted-dashed red line 
shows the correlation assuming a 0.04 fractional error on corrected distances. 


butions, mean and median values are typical for all the mocks built 
out of the different constrained simulations. Distribution of radial pe¬ 
culiar velocities can now be approximated by a Gaussian with variance 
300 km s“* at the average confidence level of 80 ± 3 % according to 
F-tests. On the other hand such a Gaussian can be rejected as a model 
at the average 98 ± 2 % confidence level for radial peculiar velocity 
distributions of biased mocks. The average F-test values are gathered 
in Table [H columns 7 and 8. 

Figure is a proof of concept that corrected distances and cor¬ 
rected radial peculiar velocities are overall closer to their true values 
than biased ones. In the top panel of this hgure, residuals between 
erroneous and true distances are drastically reduced when comparing 
corrected and true distances instead of biased and true distances. The 
standard deviation decreases from about 15 to 4 h^'Mpc. As for ve¬ 
locities, the correlation between corrected and true velocities is closer 
to the ideal one than without correction. Moreover, no additional bias 
is created since, corrected distances are not always higher/smaller than 
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original ones. Namely, distances are not all under/over-estimated. The 
reconstruction technique is put into practice in the next subsection on 
the minimized biased mocks. Fractional errors of 4% on distances are 
assumed, in agreement with the upper limit of fractional error medi¬ 
ans found in the various corrected mocks. As shown in Figure]^ this 
assumption improves the correlation between corrected and true veloc¬ 
ities. 

4.3 Reconstruction using Mocks 

For purposes of completeness, to check the robustness of the correction 
method and the independence of its outcome on the error realization, 
the halo selection and the constrained simulation, the following selec¬ 
tions are made. 

• Ten mock catalogs built out of the same distribution of halos are 
randomly selected. They differ only by the realization of errors. They 
are denoted errorMock. Namely, error distributions are still Gaussians 
with 0.17 magnitude scatter but built out of different seeds such that 
different errors are attributed to the same halo in two different mocks 
of this series ; 

• Ten mock catalogs built out of the same reference simulation are 
randomly selected. They, hereafter haloMock, differ by the selection 
of halos. Namely, halo distributions are still mimicking cosmicflows-2 
but different halos are randomly picked to build the mocks ; 

• Five mock catalogs out of five different realizations/reference sim¬ 
ulations are randomly selected. They differ by the random component 
used to build the constrained initial conditions. These five different re¬ 
alizations, hereafter referred to as simuMock, were chosen to resem¬ 
ble the Local Universe but also to present a variety of velocity field 
monopole term (term responsible for an infall if highly negative or an 
outflow if highly positive). The importance of this term will be devel¬ 
oped later in this paper. 

Next, the WF technique is applied to each one of the mocks, 
as described above, plus their corrected versions and the original 
versions. In other words, the WF is applied to; 1) the biased catalog 
(either an errorMock, a haloMock or a simuMock), 2) the corre¬ 
sponding corrected catalog (the biased mock corrected with the 
method) and, as a control, 3) the catalog with original positions and 
radial peculiar velocities, to remove any bias due to an increasing 
smoothing with distance. This sample of 66 mocks (25 biased mocks, 
25 corrected mocks and 16 original mocks, the ten mock catalogs 
built out of the same distribution of halos have the same original 
mock) are all used to reconstruct the overdensity and velocity fields 
of the reference simulations which resemble the Local Universe. 
Reconstructions obtained with the original mocks represent the best 
overdensity and velocity fields we can expect from the WF method 
for a given reference simulation and halo selection. As the goal of 
the paper is not t o test the accuracy of the WF technique, already 


widely tested (e.g. Hoffmanll 1994: Zaroubi et al.ll 19951 

1999|; 

HoffmanI 

1 19991: IZaroubiil2000l;ICourtois et al.ll2012l:ISorce et alJ 

2014a 

), it is in 


a sense easier to compare reconstructions obtained with biased and 
corrected mocks to those obtained with original mocks to determine to 
which extent biased mocks are properly corrected. For completeness 
however, the reference simulation, selected for plots in this paper, 
is shown with the three reconstructions on Figure |7] The reference 
simulation, from which the mock catalog was extracted, is plotted 
in the first column. The three reconstructions obtained with the WF 
applied to the three different mocks (original, biased and corrected) 
are shown in the next three columns. The reconstruction obtained 
with the biased mock (third column) presents what is expected from 
biases’ effects, namely a large infall in contradiction with the reference 


simulation (first column). In addition, structures are extended and very 
round. On the other hand, the strong infall has disappeared from the 
reconstruction resulting from the corrected mock (fourth column). 
Structures are more sharply defined in good agreement with the best 
result we can obtain using the WF technique on a catalog similar to 
cosmicflows-2 but original (second column), namely without errors. 
We emphasize that the WF reconstruction obtained with this original 
mock is that which reconstructs the best the reference simulation for a 
given halo selection. Other reconstructions (from mocks with errors) 
cannot be any better. 

After this qualitative analysis, a first quantitative analysis can be 
directed to the aspect of the infall. The concept is to compute monopole 
terms at 150 h^'Mpc of the WF reconstructed velocity fields, terms 
which are proportional to the divergence of these latter. Decreasing (in 
absolute value) the monopole t erm (trace of the symmetric component 
of the deformation tensor, e.g. lUibeskind et aT]| 20131 . [20141) is equiva¬ 
lent to reducing the infall onto the Local Volume. This infall is due to 
biases and as such is unphysical. The limit of 150 h^'Mpc is the result 
of the assumption that, at such a distance, the effect of radial pecu¬ 
liar velocities on galaxy total velocities has to be low. The expansion 
is predominant at such distances. As a matter of fact, monopole terms 
of velocity fields obtained with the WF technique applied to corrected 
mocks happen to be divided approximately by a factor ten in compar¬ 
ison with those of velocity fields obtained with biased mocks. Table 
[T] gathers the average monopole terms for the three different types of 
mocks as well as for the entire set of mocks. On Figure [T] the infall 
observed in the displayed reconstruction obtained with the, selected 
for display, biased mock is well reflected by the average large negative 
monopole term value (column 3: -792 ±71 km s“*). On the other hand, 
original and corrected mocks result in reconstructions without a large 
infall in agreement with low monopole term values found for the cor¬ 
responding reconstructed velocity fields and for that of the reference 
simulation (columns 2 and 4: -70 ± 47 km s“* and -58 ± 42 km s“*). 
The small scatters show that it is a general trend. Monopole term values 
are more similar between velocity fields obtained with original and cor¬ 
rected mocks than with biased mocks. The largest difference between 
monopole term values of velocity fields from reference simulations - 
or obtained with the WF applied to original mocks - and those ob¬ 
tained with corrected mocks is observed for the two simulations with 
less ‘Local Universe-like’ values (quite large absolute values, -210 and 
210 km s“'), thus they are supposed to be less representative of the 
Local Universe. As for the most common difference values, they are of 
the same order of magnitude (~ 20-30 km s“*) as the value variations 
found for velocity fields obtained with original mocks built from the 
same reference simulation but consisting of a different halo selection 
(haloMock category). 

A second quantitative analysis can be derived from cell-to-cell 
comparisons between reconstructed velocity fields obtained when ap¬ 
plying the WF technique to original and other mocks. In Table [T] 
columns 11 and 12, 1-cr scatters are smaller, by ~ 10 - 20 km s“*, 
when comparing grid cells within 320 h^'Mpc of the WF reconstructed 
velocity fields obtained with original and corrected mocks than those 
obtained with original and biased mocks. However, on Figure]^ left- 
hand panel, because the WF smooths by definition, a tilt can be ob¬ 
served in the cell-to-cell comparison of reconstructed velocity fields 
obtained with corrected and original, chosen for display, mocks. Such 
a tilt can be observed in every comparison between a velocity field ob¬ 
tained from a corrected mock and its corresponding original mock’s 
WF velocity field. Additional cell-to-cell comparisons between recon¬ 
structed divergent (due solely to densities in the box) velocity fields 
reveals that the correction combined with the WF technique affects di¬ 
vergent velocity fields, which are not affected as much as full (diver¬ 
gent plus tidal) fields by biases as shown in Table [T] columns 11 and 
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12. In other words, input velocities have to he high enough in order to 
he sufficiently high in the final reconstruction because of the Wiener- 
Filter smoothing. Currently, reconstructed velocities are too low in ab¬ 
solute value by an average factor of 1.5. A test was conducted accord¬ 
ingly. WF resulting velocities are multiplied by 1.5 in all the recon¬ 
structions obtained with corrected mocks. In every case, l-cr scatters 
are decreased some more (by ~ 5 km s“*) in cell-to-cell comparisons 
between adjusted by 1.5 velocity fields obtained with corrected mocks 
and those obtained with original mocks. Because WF reconstructed 
velocity fields abide by a prior (cosmology), this test serves the sole 
purpose of demonstrating that the tilt with an empirical 1.5 factor is 
a constant. Looking back to Figure [h] the linear fit obtained assum¬ 
ing a 0.04 fractional error on corrected distances gave us a hint about 
this 1.5 factor as, for instance, a true radial peculiar velocity of 2000 
km s“‘ correlates via the fit to a corrected radial peculiar velocity of 
1000 km s“*. Then, to remedy to the problem, corrected velocities are 
multiplied by 1.5 before the reconstruction and fractional errors are re¬ 
adjusted to 5% before applying the WF on the ’new’ corrected mocks, 
hereafter denoted 1.5-Corrected. 

The last column of Figure [7] shows the WF reconstruction ob¬ 
tained out of the randomly selected for display 1.5-Corrected mock. 
Comparing this reconstruction with the reconstruction obtained with 
the corrected mock given in the fourth column of the same Figure re¬ 
veals no major qualitative differences. Quantitatively, Figureshows 
that no evident tilt can be observed anymore when comparing the 
WF full reconstructed velocity fields obtained with original and 1.5- 
Corrected mocks. Again, in Table [T] l-cr scatters are smaller, by 20 - 
30 km s“', when comparing full velocity fields reconstructed with the 
WF applied to original and 1.5-Corrected mocks than applied to orig¬ 
inal and biased mocks. Monopole terms (in absolute value) computed 
for these reconstructed velocity fields are also decreased with respect to 
velocity field monopole terms obtained with biased mocks (by a factor 
> 5) . The average monopole terms are also gathered in Table [T] col¬ 
umn 5. The additional comparisons between reconstructed divergent 
velocity fields reveal that the correction combined with the WF tech¬ 
nique no longer affect the divergent velocity field as scatters are either 
unchanged or slightly decreased. The only (although non significant) 
exception appears for one mock of the errorMock category but this 
mock is not representative of cosmicflows-2 because of a ’relatively 
low’ monopole term (-510 km s“') in the reconstruction obtained with 
the biased mock (namely a smaller bias effect). As a matter of fact 
applying the WF to cosmicflows-2, hereafter CF2-biased, we found a 
monopole term value at 150 h^'Mpc close to -800 km s“' for the re¬ 
sulting reconstructed velocity field. 

Differential plots allow another verification which consists in 
checking that no patterns are created by the correction in WF recon¬ 
structions. These plots represent the difference or residual between 
two reconstructions. Figureshows such plots for the XY slice of the 
WF reconstructed overdensity and velocity fields obtained with the, se¬ 
lected for display, biased, corrected and original mocks. Arrow lengths 
are proportional to the residual of the velocity fields, dashed and solid 
black contours represent the residual of the overdensity fields. Green 
contours separate positive and negative residuals. From the two first 
columns, the difference between velocity field reconstructions from bi¬ 
ased/original (left) and corrected/original (middle) mocks clearly con¬ 
firms that the infall is suppressed by the correction, the residual is de¬ 
creased and no pattern is created. Numerically, the root mean squares 
of velocity field residuals (sum in quadrature of vxi — vx 2 , vyi — vy 2 , 
vzi - VZ 2 components where 1 and 2 stand for WF reconstructed veloc¬ 
ity fields obtained either with biased, original or 1.5-Corrected mocks) 
are 350 km s“' for biased/original mocks and 250 km s“* for 1.5- 
Corrected/original mocks. For comparison, the root mean square of 
the velocity field residual for two original mocks from the same refer¬ 
ence simulation (haloMock category) is about 150 km s“‘. No created 


structures are evident in the density field residuals, confirming the ab¬ 
sence of l-cr scatter changes in cell-to-cell density field comparisons, 
gathered in Table [T] The small variation in the l-cr scatters obtained 
with cell-to-cell comparisons of divergent velocity fields (derived from 
densities) constitutes another proof. On that same Figure, velocity field 
residuals, between reconstructions obtained with biased/original (first 
panel) and biased/1.5-Corrected (last panel) mocks, are similar. This 
highlights again that the correction results in reconstructions close to 
the best possible ones. These observations are valid for all the tested 
mocks. 

All these comparisons reveal that the best reconstructions (after 
reconstructions obtained with original, i.e. without error mocks) are 
obtained with 1.5-Corrected mocks. They enable reconstructions close 
to the best reconstructions possible at less than 100 km s“* (2 h^'Mpc 
in terms of displacement). For completeness, WF full velocity fields 
obtained with the original mocks are compared with velocity fields of 
reference simulations. This allows a measurement of the total effec¬ 
tiveness of the method including the WF technique applied to a mock 
similar to cosmicflows-2. Table [T] gathers the average scatter around 
the 1:1 linear relation in column 10. A sum in quadrature of the l-cr 
scatters obtained with a comparison between velocity fields obtained 
from original mocks / reference simulations and from original / 1.5- 
Corrected mocks establishes that reconstructed velocity (displacement) 
fields obtained with 1.5-Corrected mocks are good at 100-150 km s“* 
(2-3 h^'Mpc), the lowest limit that can be reached using linear theories. 


5 APPLICATION TO COSMICFLOWS-2 CATALOG 

To perfect the tests, we apply the correction, then the WF, to 
cosmicflows-2, and evaluate the monopole terms of the different re¬ 
constructed velocity fields. Figure |4] shows that the correction enables 
us to recover a close to a Gaussian radial peculiar velocity distribution 
for cosmicflows-2. Actually, a F-test confirms that a Gaussian model of 
variance 300 km s“‘ cannot be rejected at more than the 6 % confi¬ 
dence level while it can be rejected at more than the 95% confidence 
level before any correction. The CF2 line in Table[T]gives the monopole 
term values which are decreased by a factor larger than 10 for the cor¬ 
rected catalog when compared to CF2-biased. Figure[T0]shows the su- 
pergalactic XY and YZ plane of the reconstructed velocity (arrows) 
and overdensity (contours) of the Local Universe obtained with CF2- 
Biased and CF2-L5-Corrected. Galaxies i n a 5 h~‘Mpc thick s lice (red 
dots) from the 2MASS redshift catalog jHuchra et al.ll201^ are su¬ 
perimposed to the reconstructions (solely for comparison purposes). 
First the observed infall in the reconstruction obtained with the CF2- 
biased catalog is not present anymore in the reconstruction obtained 
with the corrected catalog in agreement with monopole values. Second 
structures are more sharply defined. To guide the reader’s eye, some of 
them are identified with blue names. For instance one can notice that 
the Coma region is more condensed or that the Great Wall (GW) is in 
better agreement with galaxies from the redshift survey. The same is 
also true with voids such as Bootes void which has more well-defined 
contours, contours which follow the galaxy distribution. 

Observing these changes, one can wonder whether t he Laniakea 
supercluster of galaxies discovered by ITuIIv et ^ ll2014h is affected. 
Before defining Laniakea, it can be noted that in their Daper. lTullv et alJ 
(I2OI4I) apply a Malmquist Bias correction to reduce the infall in the re¬ 
construction obtained with the cosmicflows-2 catalog. They show that 
their correction affects only the global infall and that structures are un¬ 
affected, i.e. Laniakea is not the result of the correction. Their correc¬ 
tion consists in changing only distances and keeping radial peculiar ve¬ 
locities unchanged. With this paper, we propose to change velocities to 
remove also the asymmetry bias. The Laniakea supercluster of galaxies 
is defined as a local basin of attraction and every flow which converges 
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Ref. Simu. Original Biased Corrected 1.5-Corrected 



X (h'^ Mpc) X (h‘^ Mpc) X (h‘^ Mpc) X (h'^ Mpc) X (h'' Mpc) 



-150-100 -50 0 50 100 -100 -50 0 50 100 -100 -50 0 50 100 150 -100 -50 0 50 100 150 -100 -50 0 50 100 150 

X (h-' Mpc) X (h‘^ Mpc) X (h‘^ Mpc) X (h'^ Mpc) X (h'' Mpc) 


Figure 7. XY (top), YZ (middle) and XZ (bottom) slices of the (over)density (black contours) and velocity (black arrows) fields from the reference simulation (left) 
and reconstructed with the Wiener-Filter technique applied to the original, ie. without error (middle left), biased (center), corrected (middle right) and 1.5-Con'ected 
(right) mocks. Green contours show the mean density. The net spurious infall onto the volume is clearly visible in the velocity field reconstracted from the biased 
mock and structures are round. On the opposite, structures are more sharply defined in the reconstructions obtained with the corrected and 1.5-Corrected mock and the 
strong infall has disappeared (for instance, the expulsion from the void in the top right comer of the XZ plane is recovered) in agreement with the best reconstruction 
possible, of the reference simulation (first column), with a mock similar to cosmicflows-2 but without error (second column). No major difference can be seen between 
the results obtained with the coiTected and 1.5-Corrected in agreement with the fact that they differ only in terms of the divergent part (velocities due solely to densities 
in the box) of their velocity fields. 


towards it, i.e. the boundary of the Laniakea supercluster corresponds 
to a zero velocity isosurface. Galaxies within that zone fall onto the 
local attractor while galaxies outside of that zone are going in opposite 
direction. This surface is derived with the divergent part of the velocity 
field, i.e. with densities solely in a 100 h^'Mpc radius sphere centered 
around [-47,13,-5] h^'Mpc. We have shown with cell-to-cell compar¬ 
isons that biases affect mostly the full reconstructed velocity field and 
that the divergent reconstructed velocity field is relatively unaffected 
by the biases. The correction technique proposed in this paper reduces 
drastically the average l-cr scatter in the cell-to-cell comparison of full 
velocity fields while the average l-cr scatter in the cell-to-cell compar¬ 
ison of the divergent parts of the velocity fields is mostly unchanged. 
The result is that the Laniakea supercluster of galaxies is globally un¬ 
changed. We compute the divergent velocity field within a 100 h^'Mpc 
radius sphere around [-47,13,-5] h^'Mpc and find back the Laniakea 
supercluster of galaxies as expected. 


6 SUMMARY & CONCLUSION 

The Local Universe is composed of a multitude of structures (nodes, 
sheets, filaments, voids) which need to be studied to better understand 
how gravitational motions induce the formation of such a complex 
environment. In that respect, galaxy distances and derived radial pe¬ 
culiar velocity catalogs constitute a valuable source of information to 
map the Local Universe. The Cosmicflows project develops such cat¬ 
alogs. These catalogs, however, suffer from biases which, when not 
corre cted for, re s ults in spurious velocity infall onto the Local Vol¬ 
ume. ITuIIv et ^ lIZOLlh themselves warn us that these issues must be 
addressed if the catalog is used to infer matter and velocity distribu¬ 
tions. These biases are due to the inhomogeneity of the Local Universe 
and to the fact that distances in megaparsecs are derived from distance 
moduli in magnitudes via a logarithmic function. As a result, although 
errors on distance moduli in magnitudes have a symmetric distribution, 
the distribution of uncertainties on distance estimate in megaparsecs is 
asymmetric. Consequently, the radial peculiar velocity distribution has 
a larger negative tail than a positive one. This asymmetric bias, in addi- 
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Vx,y^ (Original) Vx,y^ (Original) Vx,y^ (Original) 

Figure 8. Cell-to-cell comparisons within 320 h“*Mpc between the Wiener-Filter reconstructed full velocity fields obtained with original and con'ected mocks (from 
left to right, Corrected, Corrected but reconstmcted velocities multiplied by 1.5, 1.5-Con'ected versus original). The gradient of grey delimits the 50, 68, 80, 95 and 
99.7 % confidence zones. The tilt with respect to the 1:1 linear relation is decreased in the cell-to-cell comparisons of the Wiener-Filter reconstructed full velocity 
fields obtained with the corrected mock but reconstructed velocities multiplied by 1.5 and with the 1.5-Corrected mock versus the original. The scatter is smaller for 
the comparison between full velocity fields obtained with the 1.5-Corrected and original mocks. 


Biased - Original 1.5-Corrected - Original Biased - 1.5-Corrected 



Figure 9. XY slices of residuals between overdensity (black and dashed contours for positive and negative difference value) and velocity (black arrows) fields 
reconstructed with the Wiener-Filter technique applied to the biased (left), 1.5-CoiTected (middle) and original mocks and to the biased and 1.5-Corrected mocks 
(right). An'ow lengths are proportional to velocity residuals. Green contours separate positive and negative density residuals. The pattern of infall present in the 
differential plot between reconstructed velocity fields obtained with biased and original mocks is absent in the differential plot between velocity fields reconstructed 
with 1.5-Corrected and original mocks. No additional patterns appear. Differences in structures are similar in both cases: neither major structures, nor major voids are 
created by the correction. The right panel confirms that the reconstruction obtained with the 1.5-Con'ected mock is close to the best possible reconstruction (from the 
original mock) as the residual between velocity fields obtained with biased and 1.5-Con'ected mocks resembles that observed for velocity fields reconstiucted with 
biased and original mocks. 


tion to Malmquist effects, inserted in the Bayesian Wiener-Filter recon¬ 
struction technique, results in a spurious infall onto the Local Volume. 

This paper proposes a correclion to minimize Ihe effecl of biases 
observed in the reconstructed fields obtained with radial peculiar veloc¬ 
ity catalogs such as the second catalog of the Cosmicflows project. The 
approach is based on the radial peculiar velocity distribution rather than 
on the radial dislance distribution. Distance estimates are thus overall 
corrected retroactively after peculiar velocities are. As radial peculiar 
velocity distribulions have been shown to be Gaussians, an iteralive 
melhod is proposed to reduce spurious non-Gaussianities thus correct 
radial peculiar velocities and as a result distances. The correction is 
developed and tested on a suite of mocks. These mocks are built out of 
constrained cosmological simulations which resemble the Local Uni¬ 
verse to limil Ihe cosmic variance. They are made to mimic as closely 
as possible the error, data and uncertainty distributions of the second 
catalog of the Cosmicflows project, so as to reproduce as realistically 


as possible biases’ effects affecting radial peculiar velocity catalogs of 
our neighborhood. 

Before Ihe addilion of errors, Gaussian dislribulions are found in 
all the mock catalogs buill oul of reference simulations. The Gaus¬ 
sian distributions are similar in terms of variance and amplitude. As a 
matter of fact, F-tests revealed thal a theorelical Gaussian of variance 
300 km s“‘ is in agreement with radial peculiar velocity distributions 
of every original (i.e. without errors) mock at the average 98 ± 2 % 
confidence level. To correcl a radial peculiar velocity from a biased 
mock, both its probability to belong to this theoretical Gaussian and its 
(weighted) uncertainty are taken into account. Corrected radial peculiar 
velocilies have dislribulions in agreemenl with the theoretical Gaussian 
at the average 80 ± 3 % confidence level according to F-tests. In com¬ 
parison, such a Gaussian model for non-corrected (i.e. biased) radial 
peculiar velocity distributions can be rejected at an average confidence 
level of 98 ± 2 %. Corrected radial peculiar velocities allow the deriva¬ 
tion of new distances, required for consistency and to take care of the 
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Catalog 

1 

Ref. Simu 
km s~^ 

2 

Original 

km 

3 

Biased 
km s~^ 

4 

Corrected 
km s“^ 

5 

1.5-Corrected 
km s“^ 

6 

Original 

% 

7 

Biased 

% 

8 

Corrected 

% 

9 

1.5-Corrected 

% 

errorMock 

-90 

-80 

-761 (-780) ± 80 

-72 (-60) ± 37 

-159 (-150) ±48 

99 

1 

80 (82) ± 3 

85 (87) ± 8 

haloMock 

-90 

-71 (-70) ± 19 

-816 (-810) ±50 

-63 (-60) ± 26 

-142 (-140) ±32 

97 (98) ± 2 

1 

81 (82) ±4 

88(88)± 8 

simuMock 

-26 (-80) ± 159 

-46 (-80) ± 107 

-804 (-830) ± 79 

-18 (0) ± 58 

-60 (-70) ± 83 

97 (98) ± 2 

1 (1)± 1 

80(80)± 2 

88 (95)± 12 

All 

-77 (-90) ± 70 

-70 (-80) ± 47 

-792 (-790) ±71 

-58 (-60) ± 42 

-133 (-140) ± 62 

98 (99) ± 2 

2(1) ±2 

80 (81) ± 3 

87 (87) ± 9 

CF2 

/ 

/ 

-819 

71 

-36 

/ 

3 

77 

94 


10 11 12 13 

Ref.Simu/ Original/Biased Original/Corrected Original/1.5-Corrected 

Original Full Full / Div d Full / Div d Full / Div d 

kms“^ kms“V kms“^ unit of density km s“V kms“^ unit of density kms“V km s“^ unit of density 

h“^Mpc h“*Mpc / h“^Mpc h“*Mpc / h“*Mpc h“^Mpc / h“*Mpc 


errorMock 

84 

89 (91) ± 5/69 (70) ±3 

0.14 

79 (78) ±2/81 (80) ±2 

0.15 

71 (70) ±3/69 (69) ±2 

0.14 


1.7 

1.8 (1.9) ±0.1/1.4 (1.4) ±0.1 


1.6 (1.6) ±0.1/1.6 (1.6) ±0.1 


1.4 (1.4) ±0.1/1.4 (1.4) ±0.1 


haloMock 

86 (87) ± 2 

97 (98) ± 5/73 (73) ± 3 

0.14 

79 (80) ±3/81 (81) ±3 

0.15 

72 (72) ±3/71 (71) ±4 

0.14 


1.8 (1.8) ±0.1 

2.0 (2.0) ±0.1/1.5 (1.5) ±0.1 


1.6 (1.6) ±0.1/1.7 (1.7) ±0.1 


1.5 (1.5) ±0.1/1.4 (1.5) ±0.1 


simuMock 

85(85)± 1 

95 (94) ± 5/77 (77) ± 4 

0.14 

83 (85) ± 5/80 (79) ± 5 

0.15 

81 (86) ±7/71 (68) ±5 

0.14 


1.8 (1.8) ±0.1 

2.0 (2.0) ±0.1/1.6 (1.6) ±0.1 


1.7 (1.7) ±0.1/1.6 (1.6) ±0.1 


1.6 (1.7) ±0.1/1.5 (1.4) ±0.1 


All 

85 (84) ± 2 

93 (92) ± 6/72 (72) ± 4 

0.14 

80 (79) ±3/81 (80) ±3 

0.15 

73 (72) ±6/70 (69) ±3 

0.14 


1.7 (1.7) ±0.1 

1.9 (1.9) ±0.1/1.5 (1.5) ±0.1 


1.6 (1.6) ±0.1/1.6 (1.6) ±0.1 


1.5 (1.5) ±0.1/1.4 (1.4) ±0.1 



Table 1. mean (median) ± standard deviation of monopole terms, confidence levels obtained with F-Tests for Gaussianity and cell-to-cell l-cr scatters computed for 
the Wiener-Filter reconstructed overdensity, full and divergent velocity fields obtained with biased, con'ected and original mocks. The first column identify the mock 
series: mocks in the errorMock section differ only on the error realization (only one original mock corresponds to that series) ; mocks belonging to the haloMock 
category have a different halo selection (only one reference simulation corresponds to that series); mocks in simuMock are built from different constrained simulations. 
’Air gathers the means, medians and standard deviations over all the mocks. ’CF2’ gives values for the observational catalog, cosmicflows-2. (1) monopole term at 
150 h“*Mpc of the velocity fields of reference simulations, km s“*, (2) to (5) monopole term value at 150 h^'Mpc of reconstructed velocity fields obtained with the 
Wiener-Filter applied to original, biased, corrected and 1.5-Corrected mocks, km s“*, (6) to (9) confidence level obtained with F-Tests for Gaussianity applied to the 
original, biased, corrected and 1.5-Corrected mock catalog. A Gaussian with a variance of 300 km s^'is used for the tests. The higher the confidence level the better 
the agreement between a Gaussian with a variance of 300 km s^'and the distribution of radial peculiar velocities, % (10) l-cr scatter of cell-to-cell comparisons in 
320 h“*Mpc between full velocity fields of the reference simulations and obtained with the Wiener-Filter technique applied to original mocks, km s“* (first line) - 
h“*Mpc (second line), (11) to (13) l-cr scatter of cell-to-cell comparisons in 320 h“*Mpc between Wiener-Filter reconstructed full, divergent velocity and overdensity 
fields obtained with original and biased, corrected, 1.5-Corrected mocks, km s“* and unit of overdensity (first line) - h^'Mpc (second line). Bold figures represent 
reconstructions the closest to reconstructions obtained with original mocks. These are reconstructions obtained with the 1.5-Corrected mocks. 


Malmquist biases. The residuals of these corrected distances with re¬ 
spect to the true distances have standard deviations of approximately 4 
h^'Mpc, compared to a ~ 15 h^'Mpc before correction. 

The Wiener-Filter reconstruction technique was applied to the 
suite of mocks. Cell-to-cel comparisons between WF reconstructions 
obtained with corrected and original (in the sense that no error was 
added) mocks reveal that scatters are decreased by 10-20 km s“' 
in comparison with cell-to-cell scatters obtained for WF reconstruc¬ 
tions obtained with biased and original mocks, however, corrected pe¬ 
culiar velocities are not high enough (in absolute value) to overpower 
the Wiener-Filter smoothing. Tilts were observed in every cell-to-cell 
comparison between reconstructed velocity fields obtained with orig¬ 
inal and corrected mocks. These tilts are due to the fact that recon¬ 
structed velocities are smaller (in absolute value) than true velocities 
by a factor ~ 1.5. Tests revealed that multiplying reconstructed velocity 
fields - obtained with corrected mocks - by 1.5 removed these tilts. Still 
Wiener-Filter reconstructions abide by a prior. Accordingly, corrected 
radial peculiar velocities are multiplied by 1.5 before being inserted 
in the Wiener-Filter together with their corresponding distances. The 
Gaussian theoretical model is still valid at the average 87 ± 9 % con¬ 
fidence level according to F-tests. A 5% fractional error in agreement 
with the upper limit of the median error values found for the different 
built mock catalogs, when comparing true and corrected distances, is 
attributed to datapoint. The Wiener-Filter technique is applied to these 


new corrected mocks. Scatters obtained with cell-to-cell comparisons 
are decreased by an additional 5 km s“* and tilts are no longer evident. 

The series of built mocks highlights the independence of the 
method outcome on the added error realization, on the datapoint selec¬ 
tion and on the constrained reference simulation. With constrained sim¬ 
ulations, resembling the Local Universe, and associated mocks, mim¬ 
icking cosmicflows-2, the method robustness and accuracy have been 
tested. The resulting reconstructed overdensity and velocity fields are 
very close to those obtained when applying the Wiener-Filter to orig¬ 
inal (i.e. without error) mocks (at < 100 km s“* or 2 h^'Mpc) and 
to velocity fields of reference constrained cosmological simulations 
(at < 150 km s“' or 3 h^'Mpc). After correction, the infall is drasti¬ 
cally reduced as monopole terms are divided by 5 - 6. These monopole 
terms are related to the flow into the Local Volume, a large positive 
value or conversely a large negative value implies a large outflow, in¬ 
fall respectively. Structures are more sharply defined. Applied to the 
observational cosmicflows-2 catalog, the method confirms its ability 
to decrease the spurious infall onto the Local volume by dividing the 
monopole term of the reconstructed velocity field by a factor ~ 10. 
The corrected cosmicflows-2 catalog distribution is in agreement with 
the theoretical Gaussian at the 94 % confidence level. Structures like 
the Great Wall and Bootes Void contours are better defined. These 
changes motivate the examination of whether the Laniakea superclus¬ 
ter of galaxies defined as a lo cal basin of attract ion and flows converg¬ 
ing onto it and discovered bv lTullv et alJ (l2014h . who chose to modify 
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Figure 10. Supergalactic XY (top) and YZ (bottom) plane of the reconstructed velocity (arrows) and overdensity (contours) fields obtained with the raw or biased 
cosmicflows-2 catalog (left) and the 1.5-corrected cosmicflows-2 catalog (right). The green color stands for the mean density. Red dots are galaxies in a 5 h“*Mpc thick 
slice from the 2MASS redshift catalog. They are superimposed to the reconstructions for comparison purposes solely. A few structures are identified such as Coma, 
Perseus Pisces (PP), Bootes Void and the Great Wall (GW). For instance, the contour of Bootes Void are better defined and the Great Wall matches more galaxies from 
the redshift catalog in the reconstruction obtained with the corrected catalog. 

distances, is affected by the correction proposed in this paper. Because 
the divergent velocity field, due to densities solely in a volume, is quite 
unaffected by the biases and because the correction technique proposed 
in this paper changes only the full velocity fields (divergent velocity 
fields are quite stable), the Laniakea Supercluster of galaxies, derived 
with the divergent velocity field computed within a 100 h^'Mpc ra¬ 
dius sphere centered on [-47,13,-5] h^'Mpc, survives the correction 
proposed in this paper. 

The correction method applied to the second catalog of the 
Cosmicflows project result in a dataset of sufficient quality to build 


constrained initial conditions, based on the process described in 
ISorce et'^ l l2014ah . of similar quality from it. Simulations of the Local 
Universe constrained by cosmicflows-2 will be released in a subsequent 
paper. 
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APPENDIX 

Let us consider a particle at Lagrangian coordinates and the coordi¬ 
nates of the grid point occupied by the particle at a time t = 0 ,Xe or Eu- 
lerian coordinates. At a later time f, the particle is located on the grid at 
Xe{1) = Xi(xe) + iI/{xe, t) where iI/(xe, t) is the displa cement field from 
the initial position. The ZeFdovich approximation (IZeFdovichlll970h 
stipulates that this displacement field can be approximated by: 

iA(x£, 0 = D+(f)^/^o(x£) (6) 

where lAoCxs) is the initial displacement field. In other words, the dis¬ 
placement field behaves similarly to the perturbation density field. Its 
direction is frozen and it grows with time. 

Because the peculiar velocity field v is the time derivative of the 
distance, D{t) = a(t)xE(l), minus expansion, axE(t) = alaD(t): 

i/'(xe, 1) dD+ dt 

v(x, t) = a{t)XE{t) = aipixE, t) = ax —— x—a = d/iA(xE, t) 

D+{t) dt da 

(7) 

where / = is the growth rate, the displacement field can be 

reached through the peculiar velocity field. The ’dot’ notation corre¬ 
sponds to the time derivative. 

Today, a is Hq. Consequently, in h^'Mpc, the conversion factor 
between the velocity field and the displacement field is 100/ where f 
is the growth rate. 
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